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Abstract 

Setting up an empirical model of optical sensing to exploit the circular Bragg phenomenon displayed 
by chiral sculptured thin films (CSTFs), we considered a CSTF with and without a central twist defect 
of 7r/2 radians. The circular Bragg phenomenon of the defect-free CSTF, and the spectral hole in the 
co-polarized reflectance spectrum of the CSTF with the twist defect, were both found to be acutely 
sensitive to the refractive index of a fluid which infiltrates the void regions of the CSTF. These findings 
bode well for the deployment of CSTFs as optical sensors. 

1 Introduction 

By means of physical vapor deposition, an array of parallel helical nanowires — known as a chiral sculptured 
thin film (CSTF) — may be grown on a substrate [1, 2]. At optical wavelengths, a CSTF may be regarded 
as a unidirectionally nonhomogeneous continuum. CSTFs display orthorhombic symmetry locally whereas 
they are structurally chiral from a global perspective [3, Chap. 9], There are two key attributes which 
render CSTFs attractive for a host of applications. Firstly, CSTFs exhibit the circular Bragg phenomenon 
(just as cholesteric liquid crystals do [4]). Thus, a structurally right/left-handed CSTF of sufficient thickness 
almost completely reflects right/lcft-circularly polarized (RCP/LCP) light which is normally incident, but 
normally incident LCP/RCP light is reflected very little, when the free-space wavelength lies within the 
Bragg regime. This property has led to the use of CSTFs as circular polarization [5], spectral- hole [6], and 
Sole [7] filters, among other applications [8, 9]. Secondly, CSTFs are porous, and their multiscale porosity 
can be tailored to allow only species of certain shapes and sizes to infiltrate their void regions [10]. This 
engineered porosity, combined with the circular Bragg phenomenon, makes CSTFs attractive as platforms for 
light-emitting devices with precise control over the circular polarization state and the emission wavelengths 
[11, 12], and optical biosensors [13, 14]. 

Further possibilities for a CSTF emerge if a structural twist defect is introduced. For example, if the 
upper half of a CSTF is twisted by it/ 2 radians about the axis of nonhomogeneity relative to the lower 
half, then the co-polarized reflectance spectrum contains a spectral hole in the middle of the Bragg regime 
[15, 16]. This phenomenon may be exploited for narrow-bandpass filtering [17], as well as for optical sensing 
of fluids which infiltrate the void regions of the CSTF [18, 19]. 

We devised an empirical model yielding the sensitivity of a CSTF's optical response — as a circular 
Bragg filter — to the refractive index of a fluid which infiltrates the CSTF's void regions, with a view to 
optical-sensing applications. We considered both a defect-free CSTF and a CSTF with a central 7r/2-twist 
defect. Our model is not limited to normal incidence but also encompasses oblique incidence, and we present 
computed results for slightly off-normal incidence, a realistic situation for sensing applications. Furthermore, 
because the material that is deposited as the CSTF is not precisely the same as the bulk material that is 
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evaporated, we use an inverse homogenization procedure [20] on related but uninfiltrated columnar thin films 
(CTFs) [21] to predict the spectral shifts due to infiltration of CSTFs. 

An cxp(— iuit) time-dependence is implicit, with uj denoting the angular frequency and i = s/—l. The 
free-space wavenumber, the free-space wavelength, and the intrinsic impedance of free space are denoted 
by ko = ojy/eojM), Ao = 27r/fco, and t/q = vMoAo' respectively, with /j,q and eo being the permeability and 
permittivity of free space. Vectors are in boldface, dyadics are underlined twice, column vectors are in 
boldface and enclosed within square brackets, and matrixes are underlined twice and square- bracketed. The 
Cartesian unit vectors are identified as u x , u y , and u z . 



2 Empirical Model 

2.1 Uninfiltrated defect-free CSTF 

Let us begin the explication of our empirical model with a defect-free CSTF with vacuous void regions; 
i.e., an uninfiltrated CSTF. The z direction is taken to be the direction of nonhomogeneity. The CSTF is 
supposed to have been grown on a planar substrate through the deposition of an evaporated bulk material 
[3]. The substrate, which lies parallel to the plane z = 0, is supposed to have been rotated about the z axis at 
a uniform angular speed throughout the deposition process. The rise angle of each resulting helical nanowirc, 
relative to the xy plane, is denoted by X- The refractive index of the deposited material — assumed to be 
an isotropic dielectric material — is written as n s , which can be different from the refractive index of the 
bulk material that was evaporated [22, 23, 24]. 

Each nanowire of a CSTF can be modeled as a string of highly elongated ellipsoidal inclusions, wound 
end-to-end around the z axis to create a helix [25, 26]. The surface of each ellipsoidal inclusion is characterized 
by the shape dyadic 



u„ u„ + 7 T u T u T + 7b u b Ub, 
wherein the normal, tangential and binormal basis vectors are given as 



(1) 



u„ = -u x smx + u z cosx 
u T = u x cos x + u z sin x 



(2) 



By choosing the shape parameters 76 > 1 and 7 T 3> 1, an aciculate shape is imposed on the inclusions. For 
the numerical results presented in §3, we fixed 7 r = 15 while noting that increasing 7 T beyond 10 does not 
give rise to significant effects for slender inclusions [26]. The helical nanowires occupy only a proportion 
/ G (0, 1) of the total CSTF volume; the volume fraction of the CSTF not occupied by nanowires is 1 — /. 

At length scales much greater than the nanoscale, the CSTF's relative permittivity dyadic may be 
expressed as 

where 2f2 is the structural period and the rotation dyadics 



£ (x) = u y u y + ( u * + u z u z ) cos x + (u z u x - u x u z ) sin x 
S_ (a) = u z u, + (u x u x + u y u y ) cos a + (u y u x — u x u y ) sin a 



(4) 



The handedness parameter /;, = +1 for a structurally right-handed CSTF, and h = — 1 for a structurally 
left-handed CSTF. The reference relative permittivity dyadic has the orthorhombic form 



= £ai u„ u„ + £61 u T u T + e c i u b u 6 . 



(5) 



The nanowire rise angle x can be measured from scanning-electron-microscope imagery. In principle, the 
relative permittivity parameters {e a i, e&i; £ci} of an uninfiltrated CSTF are also measurable. However, in 
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view of the paucity of suitable experimental data on CSTFs, our empirical model relies on the measured 
experimental data on the related columnar thin films (CTFs). In order to deposit both CSTFs and CTFs, 
the vapor flux is directed at a fixed angle \ v with respect to the substrate plane. The different morphologies 
of CSTFs and CTFs are due to the rotation of the substrate for the former but not for the latter. The 
parameters {e a i, £bi, e c i, x} are functions of Xv- 

The nanoscale model parameters {n s , /, jb} are not readily determined by experimental means. However, 
the process of inverse homogenization can be employed to determine these parameters from a knowledge of 
{eai) £bi) £ci}) as was done for titanium-oxide CTFs in a predecessor paper [20]. 



2.2 Infiltrated defect-free CSTF 

With optical-sensing applications in mind, next we consider the effect of filling the void regions of a defect-free 
CSTF with a fluid of refractive index ri£. This brings about a change in the reference relativity permittivity 
dyadic. The infiltrated CSTF is characterized by the relative permittivity dyadic e , which has the same 
eigenvectors as e but different eigenvalues. Thus, the infiltrated CSTF is characterized by (3) and (5), but 
with e 2 in lieu of e , e r ^ in lieu of £^ and {e a 2, £62, ^2} in lieu of {e a i, e&i, e c i}- The nanowire rise angle 
X remains unchanged. 

In our model, the Bruggeman homogenization formalism is applied in its usual forward sense [27] to 
determine {e a 2, £b2, £c2}, from knowledge of the nanoscale model parameters {n Sl /, 7^} together with {ni, j T }, 
as described elsewhere [26]. 



2.3 Boundary- value problem 



Let us now suppose that a CSTF occupies the region < z < L, with the half-spaces z < and z > L 
being vacuous. An arbitrarily polarized plane wave is incident on the CSTF from the half-space z < 0. Its 
wavevector lies in the xz plane, making an angle 8 G [0,7r/2) relative to the +z axis. As a result, there is 
a reflected plane wave in the half-space z < and a transmitted plane wave in the half-space z > L. Thus, 
the total electric field phasor in the half-space z < may be expressed as 



E(r) 



lUy - P + lUy + P + 

a L — a R 



V2 



V2 



V2 n V2 

while that in the half-space z > L may be expressed as 



exp {inx) exp (ik$z cos ( 



exp (inx) exp (—ikaz cos 6) , z < 0, 



E(r) 



V2 



t R - 



V2 



exp (inx) exp [iko (z — L) cos 8] , z > L, 



(G) 



(7) 



wherein p± — ^fu x cos 8 + u z sin 8 and re = fco sin 8. 

Our aim is to determine the unknown amplitudes tl and r R of the LCP and RCP components of the 
reflected plane wave, and the unknown amplitudes fx and of the LCP and RCP components of the 
transmitted plane wave, from the known amplitudes and an of the LCP and RCP components of the 
incident plane wave. As is comprehensively described elsewhere [3], this is achieved by solving the 4x4 
matrix/4- vector relation 

[fexttj = \M(L)] • [f entr V]. 



(8) 



Here, the column 4- vectors 
/ 

1 

72 



(r L + rR) + {a L + a R ) 
i [- (r L - rR) + (a L - a R )] 
—i [(ri - rR) + (cil - a R )] /rj 
\ -[(r L +rR)-{a L + a B )}/vo ) 



1 

3 



t L +tR 
i {t L - tR) 
-i (t L - tR) /r]o 
(t L + tR) /770 



(9) 
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arise from the field phasors at z = and z = L, respectively. The optical response characteristics of the 
CSTF are encapsulated by the 4x4 transfer matrix \M_{L)], which is conveniently expressed as [28] 



IK(L)} = \B(h^)} • [M'(L)], 



wherein the 4x4 matrix 



cos a 


— sin a 








sin cr 


cos a 














cos a 


— sine 








sin a 


cos a 



lR(a)} = 



The 4x4 matrizant [M'(z)] satisfies the ordinary differential equation 



Tz [M!(z)} = ilg(z))>\M!(z) 



(10) 



(11) 



(12) 



subject to the boundary condition [M'(0)] = [/], with [J] being the identity 4x4 matrix. The 4x4 matrix 

-ih% 





\P'(z)] 







U)fJ,Q 



-we e c 








K,(€h— e Q ) >■ • n 

■ 2 e.T COSgSUl2x 



2e„T sing sin 2x 



- £ — sin g cos g 



- -S— sin g 



-S— cos 2 g 



sing cos g 



containing 



g = hirz/n, 



— — sin g cos g 
-2 — sin 2 g 

o 



Kief, — e a ) j- • o 

■ 2 £„t sing sin 2x 



cos 2 x + (e 6 /e ) sin 2 x 



* — cos 2 g 





(13) 



(14) 



depends on whether the CSTF is uninfiltrated or infiltrated. Equation (12) can be solved for [M'(z)] by 
numerical means, most conveniently using a piecewise uniform approximation [3] . 

Once [M'( z )] is determined, it is a straightforward matter of linear algebra to extract the reflection 
amplitudes Tl,r and transmission amplitudes th.R from (8), for specified incident amplitudes cll,r- Following 
the standard convention, wc introduce the reflection coefficients Vll,lr,rl,rr and transmission coefficients 

tLL,LR,RL,RR per 



TR 



r LL r LR 
Trl r RR 



<l fl 



tl 

tn 



tLL thR 
t R L t RR 



(I PI 



(15) 



The square magnitude of a reflection or transmission coefficient yields the corresponding reflectance or 
transmittance; i.e., R a p = \r a /3\ 2 and T a p = \t a p^ , where a,J3 € {L,R}. 



2.4 CSTF with central vr/2-twist defect 

As mentioned in §1, the introduction of a central twist defect leads to a narrowband feature that can be 
very useful for sensing applications. Therefore, we further consider the CSTF of finite thickness introduced 
in §2.3 but here with the upper half z G [L/2,L] of the CSTF twisted about the z axis by an angle of 7r/2 
radians with respect to the lower half z £ [Q,L/2). 
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Mathematically, the central 7r/2-twist defect is accommodated as follows: The relative permittivity dyadic 
of the centrally twisted uninfiltrated CSTF is per (3) but with the rotation dyadic 5 (/iff ) therein replaced 
hyg z {h[^ +V(z)]}, where 

f 0, < z < L/2 

if,(z) = I . (16) 

1 tt/2, L/2<z<L 

The relative permittivity dyadic of the centrally twisted and infiltrated CSTF follows from the corresponding 
dyadic for the centrally twisted and uninfiltrated CSTF, in exactly the same way as is the case when the 
CSTF is defect-free. The calculation of the reflectances and transmittances follows the same path as is 
described in §2.3 with the exception that (10) therein is replaced by 



{M(L)} = [B(h 



= v " L =V 2VL 2 



TlL 7T 



ttL 



[M!(L/2)}. 



(17) 



3 Numerical results 

In order to illustrate the empirical model, we chose a CSTF of thickness L = 4011 where the structural 
half-period f2 = 185 nm. The chosen relative permittivity parameters, namely 



e a i 



1.0443 + 2.7394 



1.6765 + 1.5649 



1.3586 + 2.1109 



2X„ 



2 X v 



2Xv 



1.3697 



0.7825 



1.0554 



2Xv 



2x t 



2Xv 



-i 2 



(18) 



with 

X = tan" 1 (2.8818 tan Xv) , (19) 

emerged from data measured for a CTF made by evaporating patinal® titanium oxide [21] . These relations — 
which came from measurements at Ao = 633 nm — were presumed to be constant over the range of wavelengths 
considered here. Values for the corresponding nanoscale model parameters {n s ,/, 7t}, as computed using 
the inverse Bruggeman homogenization formalism [20], are provided in Table 1 for the vapor flux angles 
Xv = 15°, 30°, and 60°. Furthermore, we set h = +1. The angle of incidence 9 was fixed at 10°. 

Computed reflectances and transmittances are plotted versus Ao in Fig. 1, for the defect-free CSTF 
for which we set Xv — 15°. Further computations (not presented here) using other values of Xv revealed 
qualitatively similar graphs of reflectances and transmittances versus Ao. The effects of three values of ni — 
namely, ng = 1, 1.3 and 1.5 — are represented in Fig. 1. The circular Bragg phenomenon is most obviously 
appreciated as a sharp local maximum in the graphs of Rrr, with attendant features occurring in the graphs 
of some other reflectances and transmittances. If A™ ax denotes the free-space wavelength corresponding to 
this local maximum, from Fig. 1, we found that \™ ax sa 622 nm for ni = 1, \™ ax f=a 712 nm for ng = 1.3, 
and Xf ax w 768 nm for n e = 1.5. 

Clearly, the circular Bragg phenomenon undergoes a substantial spectral shift as ne increases from unity. 
In order to elucidate further this matter, we focused on the spectral-shift sensitivity d\™ ax /dn,£. Graphs 
of d\™ ax /dni against A™ ax , computed for the range 1 < rig < 1.5, are presented in Fig. 2. In addition to 
results for the vapor flux angle Xv = 15°, results arc also plotted in Fig. 2 for Xv = 30° and Xv = 60°. 
For all vapor flux angles considered and all values of rig £ (1, 1.5), the spectral-shift sensitivity d\™ ax j dni 
is positive- valued and greater than 118 nm per refractive index unit (RIU). When \v = 15°, d\™ ax j dni 
generally decreases as \™ ax increases. A similar trend is exhibited for Xv = 30°, but d\™ ax / dni generally 
increases as \™ ax increases for Xv = 60°. 
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The center wavelength of the circular Bragg regime has been estimated as [29] 



A Bl ' « n + « / . ? £a2£b2 ^1 vW (20) 

V V £ a2 COS z X + £ b2 Sin X/ 

The graphs of d\$ r /dni versus Xq t , as provided in Fig. 3 for the vapor flux angles Xv = 15°, 30°, and 60°, 
are remarkably similar (but not identical) to the graphs of d\™ ax / dne, versus \™ ax displayed in Fig. 2. Thus, 
the center-wavelength formula (20) can yield a convenient estimate of the spectral-shift sensitivity, without 
having to solve the reflection-transmission problem. 

We turn now to the CSTF with a central twist defect of tt/2 radians, as described in §2.4. Graphs of 
the reflectances and transmittances versus Ao for Xv = 15° are provided in Fig. 4. As we remarked for the 
defect-free CSTF, graphs (not presented here) which are qualitatively similar to those presented in Fig. 4 
were obtained when other values of the vapor flux angle Xv were considered. The graphs of Fig. 4 are 
substantially different to those of Fig. 1: the local maximums in the graphs of Rrr in Fig. 1 have been 
replaced by sharp local minimums in Fig. 4. These local minimums — which represent an ultranarrowband 
spectral hole — arise at the free-space wavelengths A™" that are approximately the same as the corresponding 
local maximums \™ ax in Fig. 1. 

The location of the spectral hole on the Ao axis is highly sensitive to m. In a similar manner to before, 
we explore this matter by computing the spectral-shift sensitivity d\™ ln / dug at each value of ng. In Fig. 5, 
d\™ m I dni is plotted against A™ m , with the spectral-shift sensitivity computed for the range 1 < m < 1.5 
and with X v = 15°, 30°, and 60°. The plots of d\^ ln /dn e versus A^ m in Fig. 5 are both qualitatively 
and quantitatively similar to those of d\™ ax /dne versus \™ ax in Fig. 2. That is, positive- valued d\™ m /dne 
generally decreases as A™ m increases for Xv = 15° and 30°, and generally increases as A™" increases for 
Xv = 60°. A similar correspondence exists with Fig. 3. 

4 Closing remarks 

Our empirical model has demonstrated that the circular Bragg phenomenon associated with a defect-free 
CSTF, and the ultranarrowband spectral hole displayed by a CSTF with a central 7r/2-twist defect, both 
undergo substantially large spectral shifts due to infiltration by a fluid. Although, owing to lack of availability 
of experimental data, we did not consider wavelength-dispersion in the dielectric properties of the material 
used to deposit a CSTF, the promise of CSTFs — with or without a structural twist — to act as platforms for 
optical sensing was clearly highlighted. Experimental validation is planned. 
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Xv 


lb 
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n s 


15° 


2.2793 


0.3614 


3.2510 


30° 


1.8381 


0.5039 


3.0517 


60° 


1.4054 


0.6956 


2.9105 



Table 1: Nanoscale model parameters 7&, / and n s for Xv — 15°, 30°, and 60°. 
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Figure 1: Reflectances and transmittances plotted against the free-space wavelength for a defect- free 
titanium-oxide CSTF; L = 40ft, fi = 185 nm, h = +1, Xv = 15°, and 9 = 10°. The CSTF is infil- 
trated with a fluid of refractive index ng = 1.0 (blue broken-dashed curves), 1.3 (red solid curves), and 1.5 
(green dashed curves). 
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Figure 2: Spectral-shift sensitivity d\™ ax / dnt plotted against A™ aa: for ni € (1, 1.5). The vapor flux angle 
Xv = 15° (blue broken-dashed curve), 30° (green dashed curve), and 60° (red solid curve). 
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Figure 3: Spectral-shift sensitivity as estimated by dX^/dn^ plotted against A^ r for ng € (1,1.5). The 
vapor flux angle \v = 15° (blue broken-dashed curve), 30° (green dashed curve), and 60° (red solid curve). 
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Figure 4: As Fig. 1 except that the upper half z E [L/2, L] of the CSTF is twisted about the z axis by it 
radians with respect to the lower half z € [0, L/2). 
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Figure 5: Spectral-shift sensitivity d\™ m /dri£ plotted against A™ m for ri£ E (1,1.5). The vapor flux angle 
Xv = 15° (blue broken-dashed curve), 30° (green dashed curve), and 60° (red solid curve). 
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